Edge states in two-dimensional electron gas with heterogeneous spin-orbit interaction 
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We show that edge states similar to those known for topological insulators exist in two-dimensional 
electron system with one-band spectrum in the presence of heterogeneous spin-orbit interaction 
(SOI). These states appear at boundaries between regions with the SOIs of different kind or between 
the regions with the SOI and without it. Depending on the system parameters they can appear 
in an energy range lying both in the forbidden and conduction bands of bulk states. The edge 
states have chiral spin texture and carry a spin current under the equilibrium. We study also the 
size quantization of the edge states in a strip structure with two boundaries to find an unusual 
dependence of the quantization energy on the strip width. 



Spin-orbit interaction (SOI) in solids generates a 
wealth of fascinating effects which are somehow con- 
nected with the presence of boundaries and interfaces. It 
is enough to mention such effects as: the quantum spin 
Hall effect, 1,2 spin accumulation in the regime of extrinsic 
spin Hall effect, 3 spin-dependent tunneling, 4,5 edge spin 
currents^ edge spin accumulation, 7 etc. 

Presently increasing attention is paid to topological in- 
sulators 8 10 which genetically originates from the SOIs. 
Topological insulators are characterized by the presence 
of edge or surface states (correspondingly in two- or 
three-dimensional cases) lying in the energy gap of bulk 
states. An important property of topological states is 
their chiral spin texture correlated with the wave vec- 
tor, due to which these states are topologically protected 
against the scattering and carry a dissipationless spin 
current. The dispersion curve of these states connects 
the valence and conduction bands so the simplest model 
describing them uses a two-band Hamiltonian 

In this paper we show that somewhat similar states ex- 
ist also in a much more simple system where they appear 
even within one-band model. We find these states in two- 
dimensional (2D) electron systems containing a bound- 
ary between regions with the SOIs of different kind, such 
as a contact of regions with the Rashba and Dresselhaus 
SOIs. In particular, the edge states exist in a contact 
of 2D regions with the SOI and the normal 2D electron 
gas without SOI. The distinctive feature of these states is 
that the dispersion curve of their spectrum goes from the 
conduction band, circumscribes a loop in the forbidden 
band and then turns to the conduction band. 

Below we study the spectra of the edge states, the 
spin density distribution and spin-currents in structures 
RSOI/N and RSOI/DSOI, where N denotes normal 2D 
electron gas, RSOI and DSOI stand for regions with the 
Rashba and Dresselhaus SOIs. In addition, we consider 
the size quantization effect on the edge states in strip 
structures DSOI/RSOI/DSOI where the central RSOI re- 
gion is of finite width. 

The 2D system with heterogeneous SOI is described 
by the Hamiltonians 
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defined in each region indexed by i with uniform potential 
Ui and SOI strength. Here rrii is the effective mass of 
electrons, p — (p x ,P y ) is the momentum, H^ D is the 
Rashba or Dresselhaus SOI Hamiltonian, 
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a x , (jy are the Pauli matrices. We use the effective mass 
approximation which is well justified since the envelope 
wave functions are supposed to vary slowly at the lattice 
constant scaled 

To be specific, we consider the case where x axis is 
normal to the boundary and y axis is parallel to it. The 
wave functions are determined by the Schrodinger 
equation with boundary conditions 5 
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where B^ = iaidy/h 2 in the case of Rashba SOI and 
B^=ipi(j x /h 2 for Dresselhaus SOI. 
General solution in z-th region reads as 
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where k y is the tangential wave vector, k^\ is the x com- 
ponent of the wave vector, defined by a characteristic 
equation of the Hamiltonian (pQ), s stands for the spin 
index, j numbers the solutions of the characteristic equa- 
tion, \j \ i s ^ ne s pi n function. 

Generally there is a set of four wave vectors k^ s , which 
are described in detail in Ref. 5. A short resume is as 
follows. The wave vectors k^\ are complex functions of 
the energy E and the tangential momentum k y . Two of 
them correspond to the states propagating or decreasing 
to the right along the x axis, other two relate to the states 
propagating or decreasing in opposite direction. In the 
energy region E < Ui — E so , all k^ s contain both real and 
imaginary parts which describe decaying and oscillating 
states. E so is the characteristic energy of the SOI: E so = 
ma 2 /(2h 2 ) for the Rashba SOI and E so = mf3 2 /(2h 2 ) 
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for the Dresselhaus SOL When E > Ui — E so , the wave 
vectors k^ s are either purely real or imaginary depending 
on the relation between E and k y . The full analysis of the 
kj % l dependence on E and k y is available in Refs I Ml5lll6l . 

To clarify if the edge states exist near the interface x = 
one needs to find the solutions satisfying the boundary 
conditions [Eq. (j3])] and vanishing at the infinity (x — »• 
±00). Dropping the terms, which do not vanish at the 
infinity, in Eq. (j4]) we obtain a system of homogeneous 
equations. The zeros of its determinant give equations 
for the edge state spectrum. 

Below we present results of the edge-state spectrum 
calculations for several structures: contacts SOI/N, 
RSOI/DSOI and strip structure DSOI/RSOI/DSOI. 

To begin, consider the SOI/N contact. The region with 
the Rashba SOI is located at x < and N region lies at 
x > 0. The energy diagram is depicted in the insets in 
Figs [TJ [21 The potential Un in N region is lower than 
that in the SOI region: Un = —U, U so = 0. 

Of most interest is the energy interval E < —E SOl 
where the electron states in the SOI region are evanes- 
cent 5 and the wave functions are 
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( = 2m so E/h 2 , k so = m so a/h 2 , m so is the effective mass 
of electrons in the SOI region. 

In the N region the wave function is 
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where g = yk 2 — (( + u)/{i : u = 2m so U/H 2 , fi = 

m S o/mN, and tun is the effective mass in the N region. 

By matching the wave functions at x = we come to 
the following condition under which the wave function 
amplitudes are nonzero: 
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Taking into account that fc, k and g are functions of ( 
and k y , Eq. ([TP]) gives the spectrum of the edge states: 
£ = (Es(ky). This equation is rather cumbersome in the 
explicit form therefore we consider two specific cases. 

In the limiting case \i = m so /mN — > 0, the edge state 
spectrum has a simple form 
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FIG. 1: (Color online). Spectrum of the edge states at SOI/N 
boundary for m so /rriN = 0.1 and m so /mN — 0.01. Thin lines 
depict the boundaries of the bulk state continuum (shaded) 
in the SOI and N regions for U = —E so . The inset shows 
schematically the potential shape, the bulk state spectra and 
the distribution of the electron density in the edge state. 



It is seen that the edge states exist below the conduction 
band bottom in the SOI region occupying the interval 
— 1.25k g Q < C < —k 2 Q . At a given energy there are two 
pairs of the edge states with different signs of the wave 
vector k y and the group velocity. If the potential of the 
N region is higher than — 1.25E so , the edge states are 
located in the forbidden band where the bulk states are 
absent. When the effective mass ratio \i is finite, the edge 
state spectrum deviates from the limiting expression (fTTj) 
as it is shown in Fig. [1] for Un — —E so . 

At the energy higher than the conduction band bot- 
tom in the SOI region, E > —E SOl the edge states ex- 
ist against the background of the bulk states. We have 
studied these states by making specific calculations for 
this energy region. Two edge states with k y < and 
k y > are found at a given energy in the interval 
— E so < E < E cr , with E cr being the energy at which 
the edge-state spectrum intersects the bulk-state spec- 
trum boundary. 

When m so = m n , the edge state spectrum has the form 
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The validity of Eq. (fT2j) is restricted by the inequality 
ky > C + u appearing from the requirement that g to 
be real. The edge states are absent when this condition 
is violated. The edge-state spectrum has two branches 
corresponding to waves propagating in opposed direc- 
tions. They are shown in Fig. [2] for the potential step 
U = 1.5E so at the interface. The edge states occupy a 
finite energy layer and a finite interval of k y . The lower 
and upper edges of this interval are determined by the 
intersection points of the edge-state spectrum with the 
boundaries of the bulk spectra in the N and SOI regions. 
No edge states exist below the conduction band bottom 
of the N region. 
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FIG. 2: (Color online). Spectrum of the edge states at the 
SOI-N boundary for m S o/mN = 1, U = 1.5E so . 



The edge-state band bottom depends on the poten- 
tial step U at the interface. The function E^U) is easy 
to find from Eq. (flQ|) and the condition k y ) = 0, 

E b = -U+ U a ° . (13) 

It is seen that E^ = —E so at U = ^ so - With increas- 
ing C7, the edge-state bottom decreases to reach the 
minimum E^ = — 1.25E so and whereupon E^ goes up. 
Thus, the maximum depth of the edge-state bottom is 
—0.2bE so below the conduction band bottom of the SOI 
region. This conclusion is easy generalized to the arbi- 
trary mass ratio fi. 

The edge state formation can be interpreted as be- 
ing a result of the lowering of the electron energy near 
the interface because of the mutual penetration of elec- 
trons from one contacting region to another. Electrons 
penetrating from the N region into the SOI region gain 
the energy since they undergo the SOI action. On the 
contrary, the electrons of the SOI region lose the energy 
while penetrating into the N region since they do not feel 
the SOI there. If m so <C raw, the electrons penetrate 
into the SOI region much deeper than into the N region. 
Hence, the gain in the energy is larger than its loss and a 
state localized near the interface can appear with energy 
lower than the conduction band bottom. 

Another 2D system, in which we demonstrate the edge 
states in the forbidden band, is the RSOI/DSOI struc- 
ture. The edge states are studied by solving Eqs. (j3]) and 
(|4|) in the same manner as described above. In the cal- 
culations we have used the expressions for wave vectors 

k^ R J and k^ D J in the Rashba and Dresselhaus regions de- 
j-> b j-> b I — I 

rived in Ref. |5|. The results are presented in Fig. [3] for 
two potential steps at the boundary. 

The edge states are seen to exist in the forbidden band 
even if the effective masses in the contacting regions are 
equal, in contrast to the case of the SOI/N system. But 
the energy interval, where the edge states are located, 
and the general view of the spectra are quite similar to 
those of the SOI/N structure shown in Fig.[TJ The origin 
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FIG. 3: (Color online). a) Edge-state spectrum in 
RSOI/DSOI structure for U D -U R = (line 1) and U d -Ur = 
0.25E SO (line 2). Shaded regions are the bulk states. The in- 
set shows the potential shape, the bulk state spectra and the 
distribution of the electron density in the edge state, b) Spa- 
cial distribution of the spin density components (S x , S y , S z ) 
and the total spin density S for the edge state with the en- 
ergy £es = —1.2k so and momentum k y = 0.95k so near the 
RSOI/DSOI boundary in the case where Ud = Ur, a = f3, 
m R — niD. 



of the edge states can be interpreted a result of the joint 
action of the SOIs in both contacting regions. 

The electron and spin densities in the edge states are 
localized near the boundary at a distance of the or- 
der of k~J~. It is worth noting that the S x component 
of the spin density is concentrated in the RSOI region 
whereas S y component is located mainly in the DSOI re- 
gion (Fig. Et>) . The spin direction is reversed with chang- 
ing the sign of k y . Therefore the edge states are chiral 
ones. They are protected against the backscattering by 
the time-reversal symmetry. In addition, the edge states 
carry a spin current under the thermal equilibrium. 17 Un- 
der the nonequilibrium conditions when a particle current 
flows along the boundary, the spin density accumulates 
in the edge states. 

Interesting effects arise when the edge states appear 
under the condition of size quantization. We have 
studed these effects by considering the strip structures 
DSOI/RSOI/DSOI and N/RSOI/N. In this case the cal- 
culations are more cumbersome since the wave functions 
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FIG. 4: (Color online). Energy of the edge-state sub- 
band bottoms Eb as functions of the strip width w in the 
DSOI/RSOI/DSOI structure with the flat potential land- 
scape, U l ^ t — Ur — Up 9ht , a = f3, ttir — ttid- In the 
insets, the electron density distribution in the edge states is 
shown for different widths of the strip. Lines 1 and 2 (dashed) 
correspond to the lower and upper subbands. 

should be found in three regions using two boundary con- 
ditions of the form of Eq. [3l We have solved this problem 
numerically. 

The main difference of the edge-state spectra in the 
strip structures, as compared with the spectrum in the 
single RSOI/DSOI contact, is the splitting of the edge- 
state band into two subbands with different distributions 
of the electron density across the strip. This is demon- 
strated in Fig. 2] where the energies E^i and of the 
edge-state band bottoms in the RSOI/DSOI/RSOI struc- 
ture are drawn as functions of the strip width w. 

It is interesting that the quantization energy depends 
on the width in unusual manner. The lower subband 
bottom Em decreases with increasing w until wk so < 1.5. 
This is trivially explained by the decrease of the kinetic 
energy. However, E^i unexpectedly grows as wk so > 1.5. 

To interpret such a behavior of the quantized energy 
let us take into consideration the fact that the electron 
density is redistributed across the strip with increasing 
w. The electron density in the lower subband is redis- 



tributed from the center to the edges. As it has been 
discussed above, electrons gain the energy near the in- 
terface because of the mutual action of the SOI in both 
regions. In the case of the strip, there is an additional 
energy gain caused by the joint effect of two interfaces. 
This energy gain decreases with w because the edge states 
overlap less. It is that reason due to which E^i grows with 
w. In the case of the upper subband this effect is much 
smaller since the electron density in these states is always 
concentrated closer to the edges. 

It is worth noting that this effect results in an essential 
lowering of the edge state energy down to the forbidden 
gap as compared with the case of a single interface. 

Above we considered the structures with step-like 
change in the SOI parameters at the boundaries. So 
the question arises whether the obtained results persist 
in the case where the SOI strength changes smoothly in 
the transition layer between the regions with different 
SOIs. We studied the problem and found that the edge 
states survive if the transition layer width L is small, 
Lk so < 1. With increasing the layer width the energy 
interval, where the edge states are located, diminishes 
to zero and the states disappear. In this point, these 
edge states differ from the edge states in 2D topological 
insulators. 

This is not only difference between these systems. In 
contrast to the topological edge states, the states stud- 
ied here exist in 2D electron gas with one-band spectrum 
and appear in a finite range of the parameters (the SOI 
strength, potential step hight, the interface structure pa- 
rameters). In addition, their spectrum has two pair of 
branches due to which the interbranch scattering is pos- 
sible. 
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It is known that there is another mechanism of the edge 

spin currents in 2D systems with heterogeneous SOI. It 

is caused by the spin-dependent scattering of bulk elec- 



trons on the interface.— These currents also exist under 
the equlibrium. The estimations show that the edge- state 
spin current is small compared to scattering spin current 
if the Fermi energy lies higher than E so above the con- 
duction band bottom. In opposed case, the edge-state spin 
current is prevailing, especially when the Fermi level lies 
in the forbidden band. 



